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Abstract. We provide a detailed quantum treatment of the spectral characteristics and of the dynamics of 
nondispersive two-electron wave packets along the periodically driven, collinear frozen planet configuration 
of helium. These highly correlated, long-lived wave packets arise as a quantum manifestation of regular 
islands in a mixed classical phase space, which are induced by nonlinear resonances between the external 
driving and the unperturbed dynamics of the frozen-planet configuration. Particular emphasis is given to 
the dependence of the ionization rates of the wave packet states on the driving field parameters and on 
the quantum mechanical phase space resolution, preceded by a comparison of ID and 3D life times of 
the unperturbed frozen planet. Furthermore, we study the effect of a superimposed static electric field 
component, which, on the grounds of classical considerations, is expected to stabilize the real 3D dynamics 
against large (and possibly ionizing) deviations from collinearity. 

PACS. 32.80.Rm Multiphoton ionization and excitation to highly excited states (e.g., Rydberg states) - 
05.45.Mt Semiclassical chaos / quantum chaos - 31.25.Jf Electron correlation calculations for atoms and 
ions: excited states 

1 Introduction a quantum wavefunction and a classical point particle j|, 

and has led to various applications in the context of quan- 
As documented by an increasing amount of literature, turn control (e.g., j|,|j|). For a long time, it was common 
electronic wave packets represent an important object of belief that the maximum evolution time until which such 
experimental studies in atomic and molecular physics (e.g., a localized phase space distribution can be maintained is 
The possibility to localize the electronic population ultimately limited by the time scale corresponding to the 
in phase space can be regarded as the crucial link between 
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anharmonicity of the underlying classical potential. Af- nomenon is best understood from the underlying classical 

ter this time scale, any coherent superposition which was dynamics. In classical phase space, the time-periodic per- 

initially well localized in position and momentum space turbation of the microwave destroys the global integrabil- 

would have spread out over the whole classical orbit due ity of the Coulomb dynamics; local "islands" of regular 

to dispersion, i.e., due to the fact that different compo- motion, embedded into a "sea" of chaotic, ionizing dy- 

nents of the superposition contribute with different clas- namics, are created by the nonlinear resonance between 

sical frequencies. This is usually termed as the collapse of the driving field and the unperturbed Kepler motion. The 

the wave packet. nondispersive quantum wave packets are nothing but the 



time-periodic eigenstates (Floquet states) of the driven 
atom, which arise from a local quantization within these 



resonance islands Q,[15| (see |16 for a recent review on 
the topic). 



A decade ago, however, it was found that this spread- 
ing can actually be suppressed by means which are com- 
mon in quantum control, namely by applying electromag- 
netic fields to the quantum configuration. This was explic- 

.,, , , r , l • i The above interpretation in terms of nonlinear reso- 

ltiy demonstrated tor one-electron atoms driven by mi- ^ 

c , i • 1 1 i • iHRiai i rRl imi j i nances in a mixed phase space suggests that the concept of 

crowave fields with linear |6|,[7],|§|, circular ||9|,|10||, and el- 00 

,. , HT-n, mi n ji i , • • nondispersive wave packets can be generalized to atomic 

hptic polarization |llll| . I here, a Rydberg wave packet, mi- ft a 

i • ii i ill tv i , ., .,, . . , (or molecular) systems that are more complex than hy- 

tiaily launched along a Kepler orbit with principal quan- v ' J r J 

, . , , c ,i ,. , drogen or hydrogen-like Rydberg atoms - such as, e.g., 

turn number n, is exposed to a microwave held which ° ° J ° ° 

.., ., . , . , .. t> the correlated dynamics in doubly excited helium. How- 

is resonant with the periodic classical motion, ihcn, for J J 

. , „ , n ,., , / i. i , ever, such a generalization is not at all obvious, since most 

appropriate held amplitudes (which must neither be too 

, , , , , ^ , , r doubly excited two-electron orbits of the bare three-body 

strong nor too weak, compared to the Coulomb force on J J 

i s , e ii • j_ i i i Coulomb problem (as the classical version of the doubly 

the electron) and tor the appropriate relative phase be- v 

. n , , i.i -rr n .,, , . excited helium atom) are unstable and promote rapid au- 

tween the microwave field and tfie Kepler oscillation, the ' 

,. r , i , , / , . , ,, toionization of the configuration. This prevents the con- 

spreadmg of the wave packet (which usually occurs after ° * 

, r „ , i ii, „ , , . , finement of an initially localized wave packet on longer 

approx. n/[6TT{An) \ Kepler cycles, given a coherent su- J ° 

perposition of An unperturbed eigenstates centered around sca l es - 

n) is completely inhibited. While perfectly keeping their Nonetheless, extensive studies on the classical dynam- 

shape, such nondispersive wave packets may well follow ics of helium JlTf have identified local regions of regular 

the classical motion over time scales up to 10 6 Kepler cy- motion which are associated with configurations of partic- 

cles ^,|||l^,|T^| before appreciable ionization (induced by ular symmetry of the two-electron atom. The most promi- 

the driving field) sets in This stabilization phe- nent and surprising example is the "frozen planet" orbit 
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|l8|,[l9|, a highly polarized, nearly collinear configuration sical considerations, at experimentally accessible quantum 
with asymmetric excitations of the electrons. The frozen numbers of doubly excited helium j2l],^|. Yet, this con- 
planet's stability essentially relies on the non-negligible jecture remains to be corroborated by an exact quantum 
electron-electron interaction (the electrons are, contrary treatment. 

to intuition, located on the same side of the nucleus). Sta- ., , , . .. , r 

' ' With the presently available computational facilities, 

ble eigenmodes sustain radial, vibration-like oscillations of r nl nn , . ... , , ,. r v. • i_i 

° tuli blown quantum ao initio calculations ot highly cor- 

the outer electron, as well as coupled transverse excursions , , , , , , , , , , , 

1 related two-electron wave packet states under clcctromag- 

of the outer and inner electron. Well-defined eigenfrequen- . , . . , TT , . 

° netic driving are prohibitive. However, since the stabilizing 

cies are associated with these eigenmodes, what makes the , c c , , , • • ., 

b ' static field confines the wave packet to the near vicinity 

configuration a good candidate for inducing nonspreading c ,, c , , , . ,. . ,. , ,. e 

° ° ° * ° of the held polarization axis, essential properties of these 

wave packets in the correlated two-electron dynamics, by . , .. . , ...... 

J J wave packets appear accessible through the restriction ot 

near-resonant external driving. , , . . T , , n . 

° the quantum dynamics along this axis. Indeed, nondisper- 

T r , . i • r ^1 i • i i c sive two-electron wave packets associated with the classi- 
m tact, our earlier analysis ot the classical dynamics oi 

, ,. . ,. , i ji 'ciii j cal resonances are found in quantum calculations on the 
helium m a linearly polarized electromagnetic held already 

, , c , , collinear model atom, as we have already reported in a 

provided strong support tor nondispersive wave packets J 

, c , FTTlCTni rrni . i previous letter p4| . Here we shall give detailed account of 

along the trozen-planet orbit |20[|21||. ihe external per- 1 i— <> ° 



, , . . n , . , n ... . , our quantum treatment of these obiects, identify them as 

turbation induces regular islands within the phase space 

r it i i i j i v j special eigenstates in the Floquet spectrum of the driven 

ot collinear motion (where electrons and nucleus are lined r ° i r 

„ , , i • c u i • \ i-i atomic system (alike their counterparts in driven hydrogen 

up parallel to the driving held polarization axis), which j \ t- jo 

, , , , . , atoms), and investigate their characteristic properties, 

are due to nonlinear resonances between the drive and " ° 

the oscillation of the outer electron. The dynamics within To make the present contribution self-contained, we 

these resonances is in general unstable with respect to start out in Section |^ with a review of the classical dy- 

deviations from collinearity (in contrast to the analogous namics of the driven frozen-planet configuration, which 

scenario in atomic hydrogen, the resonantly driven one was already described elsewhere p0| , [2l| ] . Section |^ is de- 

dimensional Kepler orbit [p2| ). However, a static electric voted to the formulation of the quantum eigenvalue prob- 

field applied parallel to the microwave polarization axis lem for the collinear model. As a first application thereof, 

can be used to enforce stability of the driven configura- we study the unperturbed collinear atom in Section ||, 

tion in all three spatial dimensions pO| . This suggests that what provides a remarkable example of regular dynamics 

nondispersive two-electron wave packets can be built on in an atomic two-electron system. In Section |5|, we show 

these stabilized resonances, and, on the basis of semiclas- how the nondispersive wave packets emerge in the spec- 
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trum of driven helium (Section 5.1) and discuss their life- formed according to 
times (Section p.2\), as well as their properties in presence 

Ti^N 2 ri (4 = 1,2), (3) 

of an additional static electric field (Section |5.3| ). 

Pi>— >iV -1 Pj (* = 1,2), 
1 1 — > N 3 t, 

2 Classical dynamics of the driven F 1 — y N~ F, 



frozen-planet configuration 



N' 3 lu, 



H i — > N H, 



2.1 The frozen-planet configuration of helium where N re p r esents an arbitrary, real positive quantity. 

Since any classical action scales linearly with N (as does 
In atomic units, which are used throughout this paper, the product of position and momentum), we shall iden- 
the classical Hamiltonian of the electromagnetically driven tify N with the principal quantum number of the inner 
helium atom reads electron, which corresponds to the action variable of its 

Kepler orbit. Due to the above scale invariance, we can 

H = pi + P| _ z_ Z_ _J_ 

2 2 | r"x | I r2 1 restrict our classical phase space analysis to a fixed quan- 

+ F(t) (ri +T2), (1) turn number N = 1, and use (||) to deduce the actual 

phase space structure of the energy range of interest. 

with Vi = (xi,yi,Zi) and Pi = {p lx ,Pi y ,Piz) the position 

Figure fy shows a solution of the classical equations 

and the momentum of electron i — 1,2, respectively, Z — 2 

of motion generated by (|]), which illustrates the frozen 

the nuclear charge, and 

planet configuration we shall focus on throughout this pa- 

„ _ . per. Both electrons are located on the same side of the nu- 

F = (F cos tot + F st )e z (2) K 

cleus, librating with respect to a common symmetry axis, 

the external, linearly polarized driving at frequency to While the inner electron oscillates on extremely eccentric 

and amplitude F, possibly superimposed by an additional, Kepler trajectories around the nucleus, the outer electron 

static component with field strength F st (e z represents is dynamically stabilized due to the rapid oscillation of 

the unit vector along the z axis). In analogy to driven the inner, and remains nearly "frozen" in the vicinity of 

hydrogen f2q] as well as to the unperturbed helium atom an equilibrium distance where the attractive and repul- 

p6fl . the Hamiltonian (|l|) exhibits general scaling laws: sive forces due to the other charged particles cancel each 

The classical dynamics generated by (Q) remain invariant other when averaged over one Kepler cycle of the inner 

if all variables and parameters of the system are trans- electron. The configuration is classically stable against au- 
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Fig. 1. A typical (planar) trajectory of the frozen-planet configuration. While the inner electron oscillates on eccentric Kepler 
trajectories around the nucleus, the outer electron is dynamically stabilized and remains nearly "frozen" around a given equi- 
librium distance. The configuration is stable against autoionization and represents a large regular region in the classical phase 
space of helium. 



toionization and defines a relatively large region of regular 
motion in the dominantly chaotic phase space of helium. 
Exact quantum calculations on doubly excited helium Jr9| 
have indeed revealed the existence of long-lived autoioniz- 
ing states (the energetically lowest one, below the N = 3 
threshold) which are localized along the frozen-planet or- 
bit. 

For a suitable choice of initial conditions, the motion 
of the electrons remains confined to the z axis, what de- 
fines the subspace of collinear motion. In this simplest 
case, the regular frozen-planet dynamics is characterized 
by two modes which almost perfectly separate (Fig. |J): 
the fast Kepler mode of the inner electron (which is also 
manifest in the motion of the outer electron, see Fig. 0(b)) 
and the slow oscillation of the outer electron around the 
equilibrium distance (which also affects the motion of the 
inner electron, as a slow modulation of its maximal excur- 
sion). 



This separation of time scales (the Kepler oscillation 
is almost 15 times faster than the slow mode of the outer 
electron) allows to treat the frozen-planet dynamics within 
the formalism of adiabatic invariants p7| . This defines an 
effective potential which describes the slow dynamics of 
the outer electron, in the combined field of the nucleus 
and of the rapidly oscillating inner electron pq] . As we 
see in Fig. ||, this potential is attractive Coulombic for 
large distances (due to the attraction by the nucleus), 
strongly repulsive for short distances (due to encounter 
with the inner electron) and exhibits a local minimum at 
the equlibrium distance. 

From the shape of this potential, we can immediately 
infer intrinsic scales for the frequency and the field strength, 
which depend only on the inner electron's quantum num- 
ber N, and which completely determine the effect of an 
external drive on the configuration. The natural scale for 
the field strength, Fj ~ 0.03 N~ A , is given by the maxi- 
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Fig. 2. Collinear trajectory of the frozen-planet configuration. The positions of the electrons are plotted as a function of time. 
We see that the dynamics is characterized by two almost separable modes: (i) the Kepler mode of the inner electron, which also 
manifests in the motion of the outer electron, as shown in the magnification (b) of the latter's trajectory; (ii) the slow oscillation 
of the outer electron around the equilibrium distance, which induces a feeble (and, at this scale, hardly visible) modulation of 
the amplitude of the inner electron's excursion. 




-0.2 



z 1 /N' 



Fig. 3. Effective adiabatic potential describing the slow oscillation of the outer electron around its equilibrium position (see 
[psj). Intrinsic scales for frequency and field strength can be infered from the shape of the potential. The natural scale of the field 
strength, Fi ~ 0.03 JV -4 , is given by the maximum slope of the potential (at z\ ~ 3.7 iV 2 ). The frequency scale, uii ~ 0.3 iV~ 3 , 
is given by the curvature of the potential at its minimum at Z\ ~ 2.6 N 2 , i.e., by the frequency of small oscillations around the 
equilibrium. 

mum slope of the potential (at Z\ ~ 3.7 7V 2 ; throughout frequency of small oscillations around the equilibrium po- 
the paper, subscript 1 denotes the outer, and subscript sition. In the following, we shall study the effect of an elec- 
2 the inner electron), and indicates the maximum static tromagnetic field which resonantly drives the outer elec- 
field that can be applied to the configuration without ion- tron's slow oscillation, i.e., a field with amplitude F < Fj 
izing it. The frequency scale u>j ~ 0.3 N~ 3 is given by and frequency u> ~ luj. 
the curvature of the potential at its minimum, i.e., by the 
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2.2 Phase space of the collinear, driven configuration cot — <f>o (mod 2ir). As a consequence, we effectively plot 

the outer electron's position and momentum for z-i = 



and for fixed tot. The reduction of the phase space to a 
two-dimensional surface is completed by restricting the 
initial conditions to those with fixed value N = 1 of the 
inner electron's action. This latter quantity represents an 
adiabatic invariant of the system and therefore remains 
nearly constant as time evolves. 



In the presence of the time-periodic perturbation, the dy- 
namics of the collinear, driven frozen-planet evolves in a 
five-dimensional phase space spanned by the positions and 
momenta of the electrons, and by the phase cut of the driv- 
ing field. A complete visualization of the driven dynamics 
within a simple two-dimensional Poincare surface of sec- 
tion is therefore not possible. However, for oj ~ uji, and for 

field amplitudes F < Fj which do not appreciably affect Fig. |J shows the Poincare section that is obtained by 

the Kepler motion of the inner electron, the separation of the two-step method described above, for fixed frequency 

time scales allows to map the phase space structure onto a u> = 0.2 iV~ 3 < uji, fixed driving phase cut — 0, and 

two-dimensional surface by means of a two-step Poincare for variable field amplitude F. In the unperturbed case 

section method [^l[|2^]. We begin with a simple Poincare (F = 0, Fig. |](a)), we recognize a regular phase space 

section by plotting the outer electron's phase space vari- structure, with closed curves corresponding to the regu- 

ables Z\, p z \ whenever the inner electron reaches the nu- lar oscillation of the outer electron within the effective 

cleus (22 = 0) - i.e., whenever the phase of the Kepler potential. If we switch on the driving field, phase space 

oscillation assumes a fixed value. In the time domain, the turns mixed regular-chaotic. The bounded frozen-planet 

points obtained by this section are separated by the Kepler dynamics is now represented by a local regular region cen- 

period which is much shorter than any other time scale in tered around the equilibrium point, outside which the dy- 

the system. The dynamics generated by the slow modes - namics is chaotic and leads to ionization. This remaining 

i.e., the outer electron's oscillation within the effective po- regular region is, due to the Kolmogorov-Arnold-Moser 

tential, as well as the external driving - is therefore very (KAM) theorem J27J, still large for weak field amplitudes 

well monitored by this sequence of points. (F = 0.001 A^ -4 , Fig. f|(b)), but decreases in size with 

Neglecting the effect of resonances between fast and increasing F. As a consequence of the Poincare-Birkhoff 

slow dynamics, the intersection of the Lagrangian mani- theorem (2?]], elliptic substructures, induced by nonlinear 

fold with the surface defined by z^ — can now be very resonances between the external driving and the unper- 

well approximated by a cubic interpolation between sub- turbed oscillation of the outer electron, emerge within 

sequent points in the above sequence. This yields a con- the regular region. The prominent substructure located 

tinuous trajectory which can be used to perform a second around z\ ~ 4.8 A^ 2 for F = 0.001 A~ 4 (Fig. |4|(b)) corre- 

Poincare section, by fixing the phase of the driving field sponds to the 1:1 resonance where one oscillation cycle of 
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Fig. 4. Phase space of the collinear, driven frozen-planet configuration, for variable field amplitudes, (a) F — 0, (b) 0.001 N~ 4 , 
(c) 0.005 iV -4 , at fixed driving frequency uj = 0.2 N~ s . The momentum p z \ of the outer electron is plotted versus its position 
z\ at fixed action N = 1 of the Kepler mode, fixed position Z2 — of the inner electron, and fixed phase wt = (mod 2ir) 
of the driving field |^l],^|. The regular phase space of the unperturbed atom (a) turns mixed regular-chaotic in presence of 
the external perturbation (b, c). The nonlinear 1:1 resonance between the driving field and the unperturbed oscillation of the 
electron gives rise to a substructure at F = 0.001 TV -4 (b) (centered around Z\ ~ 4.8 N 2 ), which transforms into a separate 
regular island within the chaotic sea at larger field amplitude F = 0.005 N~ 4 (c). 



the outer electron is completed after precisely one period 
of the driving field. 

At larger field amplitude (F ~ 0.003 A^ 4 for the fre- 
quency used in Fig. ||), the field-induced resonance even- 
tually disconnects from the rest of the regular domain. 
The phase space then exhibits two separate, large regu- 
lar islands embedded into the chaotic sea (see Fig. ||(c) 
at F — 0.005 A -4 ): the intrinsic island (centered around 
Z\ ~ 2.5 N 2 in Fig. [|(c)) which is essentially due to the 
intrinsic nonlinear frozen-planet dynamics with the ex- 
ternal field acting only as a small perturbation, and the 
field-induced 1:1 resonance island (centered around z\ ~ 
4.9 A" 2 in Fig. 0(c)) which arises from the combination of 
two equally important nonlinear components: the coupling 
to the external driving field and the internal interactions 
between the charged particles of the atom. 



Fig. U illustrates the effect of tuning the frequency 
at fixed field amplitude F = 0.005 A^ 4 . With increas- 
ing w, the field-induced resonance shifts towards lower ex- 
citations within the effective potential, i.e., towards un- 
perturbed orbits with higher frequencies. Consequently, 
the resonance island is located closer (as compared to 
Fig. H(c)) to the equilibrium distance for uj = 0.25 A -3 
(Fig. ^|(a)) and finally undergoes a smooth transition into 



the intrinsic island at to — 0.3 A^ 



(Fig. |(b)). 



Above the intrinsic frequency loi - which corresponds to 
the maximum frequency of free oscillations within the ef- 
fective potential - the electromagnetic field can no longer 
induce a 1:1 resonance. Resonances of higher order, how- 
ever, are still possible (e.g., the 2:1 resonance, where one 
period of the outer electron is completed precisely after 
two field cycles; the elliptic islands created by this reso- 
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Fig. 5. Phase space of the collinear, driven frozen planet configuration for variable driving frequency, (a) lo — 0.25 N 3 , (b) 
0.3JV" 3 , (c) 0.4 iV~ 3 , at fixed amplitude F = 0.005 N~ 4 and phase ut = of the drive. With increasing lo, the 1:1 resonance 
island (centered around Z\ ~ 3.9 7V 2 in (a)) shifts towards the equilibrium distance, and finally undergoes a smooth transition 
into the intrinsic island at lo — uii — 0.3 TV -3 (b). The two islands at z\ ~ 4 N 2 , p z \ ~ ±0.2 TV -1 for lo = 0.4 iV~ 3 (c) are induced 
by the 2:1 resonance where one period of the outer electron's motion matches two field cycles. 



0.6 




-06 1 ' 1 ' - ■ '■• f r ' 1 1 ' • 1 I • ' I 

5 5 5 
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Fig. 6. Phase space of the collinear, driven configuration for variable phases of the driving field, (a) cut = 0, (b) 7r/2, (c) it, 
at fixed frequency lo — 0.2 N~ 3 and field amplitude F — 0.005. While the intrinsic island remains basically at rest during time 
evolution, the field-induced 1:1 resonance island faithfully traces the resonantly driven trajectory of the outer electron. 

nance are located at z\ ~ 4 N 2 and p z \ — ±0.2 iV -1 in during the time evolution. The field-induced 1:1 resonance 
Fig. H(c)). island, on the other hand, oscillates once around the in- 

trinsic island within one field cycle. 

The crucial qualitative difference between the intrinsic 
island and the field-induced resonances is best illustrated 

2.3 Stabilization against deviations from collinearity 

by visualizing the phase space structure for different times, 

i.e., by performing the above two-step Poincare section for In contrast to their counterparts in driven hydrogen [ fi6[ 
different phases cot of the driving field. Fig. |] shows the ^9|, however, the resonance islands of the collinear, driven 
phase space of the driven configuration at lo — 0.2 N~ 3 , frozen-planet dynamics are not embedded in a regular 
F = 0.005 TV -4 , for the field phases tot — 0, 7r/2, and tt. phase space domain of the three-dimensional problem. 
We see that the intrinsic island remains basically at rest Fig. |?| shows the evolution of the configuration launched at 
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Fig. 7. Unstable trajectory of the driven frozen-planet configuration in two-dimensional configuration space. The configuration 
is initially launched at the center of the collinear 1:1 resonance island for uj = 0.25 iV" 3 and F = 0.005 N~ 4 (see Fig. |(a)), 
except for a small transverse component yi = 0.01 N 2 of the outer electron's position. The resulting time evolution of the 
trajectory is depicted from to 15 field cycles in (a), from 15 to 30 field cycles in (b), and from 75 to 91 field cycles in (c). We 
see that the initially small deviation of the configuration from the z axis increases with time, until, after approx. 10 field cycles, 
both electrons flip over to the other side of the nucleus (a). The configuration then performs several chaotic rotations around 
the nucleus (b), before the correlation between the electrons eventually breaks down and the atom ionizes (here at about 90 
field cycles) (c). 



the center of the 1:1 resonance island of Fig. |](a) (i.e., for 
uj = 0.25 N~ 3 and F = 0.005 iV" 4 ) if the outer electron 
is slightly displaced from collinearity by a small trans- 
verse component y-y — 0.01 N 2 of its position. Due to the 
effective torque exerted by the driving field on the config- 
uration, the initially small deviation from the field polar- 
ization axis increases rapidly with time, until, after about 
ten field cycles, both electrons "flip over" to the other 
side of the nucleus. The configuration then performs sev- 
eral chaotic rotations around the center, until the correla- 
tion between the electrons eventually breaks down and the 
atom ionizes after a subsequent electron-electron collision. 
This transverse instability turns out to be characteristic 
not only for the 1:1 resonance, but also for the intrinsic 



island, as well as for the 2:1 and for higher resonances (ex- 
ceptions are tiny resonances of higher order, such as the 
3:2 resonance, see |2l|,|^|). To launch nondispersive wave 
packets along the frozen-planet orbit of the real three- 
dimensional atom, an additional perturbation of the atom 
is therefore required, in order to stabilize the electronic 
motion with respect to deviations from collinearity. 

In fact, such a stabilization can be achieved by adding 
a static component F st to the electromagnetic field (see 
Eq. @). In resonantly driven hydrogen, a static field, ap- 
plied parallel to the microwave polarization axis and forc- 
ing the electron away from the nucleus, tends to confine 
the electron to the immediate vicinity of the polarization 
axis This confinement effect can now be im- 
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Fig. 8. Stabilization of the transverse frozen planet dynamics at uo — 0.25 TV -3 and F = 0.005 iV -4 . As in Fig. |^, the configuration 
is initially launched at the center of the collinear 1:1 resonance island, except for a small transverse component j/i = 0.01 iV 2 
of the outer electron's position. Without the static field (a), the electrons drift away from the z axis and enter into chaotic 
rotations around the nucleus. A static electric field, F st = 0.1F, along the z axis, directed such as to force the electrons into 
the positive z direction, counterbalances this drift mechanism (b); the electrons are driven back towards the z axis, and the 
configuration performs stable, quasi-periodic oscillations around the fundamental orbit shown in (c). 



ported to stabilize the motion of the driven frozen-planet 
- by preventing the configuration from chaotic rotations. 
Indeed, for weak static fields, F st ~ 0.1F, and for not 
too strong driving amplitudes F (e.g., F = 0.005 iV -4 at 
uj = 0.25 7V~ 3 , see Fig. small deviations from collinear- 
ity do not lead to ionization. Instead of drifting away to 
the other side of the nucleus (as in the absence of the static 
field, Fig. |§](a)), the electrons are forced back towards 
the field polarization axis (Fig. ||(b)), and the configu- 
ration performs stable, quasi-periodic oscillations around 
the fundamental orbit shown in Fig. §(c). 

This stabilization mechanism permits to embed the 
nonlinear resonances into regular regions within the un- 
restricted phase space of three-dimensional motion, while 
the electronic dynamics within the collinear subspace is 
only marginally modified by the weak static field. We find 
that this multidimensional phase space region is partic- 



ularly large for the 2:1 resonance (Fig. ||(c)): By means 
of the semiclassical Einstein-Brillouin-Keller (EBK) quan- 
tization criterion [fn) (which typically yields rather pes- 
simistic results, see |3^]), we estimate that the minimum 
principal quantum number of the inner electron needed to 
fully localize a quantum state on the stabilized 2:1 reso- 
nance island lies in the range N ~ 50 . . . 100 pl| , p3| . Such 
double excitations are clearly beyond reach in state-of- 
the-art experiments on helium, which currently can ex- 
cite N ~ 10 [ p3| . However, they should be accessible in 
earth alkaline atoms such as barium, where multi-step 
isolated-core excitations |34|] permit selective transitions 
up to states with quantum numbers N ~ 100 of the inner 
electron |3j|. At N — 50, the field parameters that are 
required to generate the 2:1 resonance islands of Fig. |^(c) 
are given by cu/(2n) ~ 20 GHz and F ~ 8 V/cm. Hence, 
by virtue of this semiclassical argument, nondispersive 
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two-electron wave packets appear realizable in earth alka- lar momenta coupled by the field, at comparable or even 
lines, at microwave frequencies and field amplitudes com- higher excitations, 
monly used in state-of-the-art ionization experiments on 

Nevertheless, essential qualitative information on the 

one-electron atoms ]36| . 

wave packets can be obtained by an approximate quan- 
tum description, which is motivated by the classical mo- 

3 Formulation of the quantum eigenvalue ^ on °^ ^ ne electrons. As discussed in the previous section, 

the static field confines the electrons to the vicinity of the 

problem 

field polarization axis (Fig. ||), without major modifica- 

_ _. ... .. „ . tions of the dynamics in the collinear subspacc. Hence, also 

3.1 The collinear Hamiltonian J r ' 

a quantum state associated with the transversly stabilized 

The preceding classical analysis strongly suggests the exis- wave packet orbit should exhibit strong localization close 

tenceof nondispersive wave packets which propagate along to the z axis. We therefore expect that essential proper- 

the frozen-planet orbit, together with a rough estimate of ties of such a quantum state are reproduced by a collinear 

the minimum required quantum number for their unam- model which restricts the electrons to move along the field 

biguous identification in the spectrum of doubly excited polarization axis. This, however, requires that the model 

helium. Reliable information on the quantum properties takes into account the full Coulomb interaction between 

of the field-induced resonances, however, can only be ob- the charged particles (i.e., that it really represents the 

tained by quantum calculations. In particular, classical exact quantum analog of the classical, collinear configura- 

mechanics tells nothing about the lifetimes of the wave tion). A smoothing of the Coulomb singularity, which is 

packets, since these are associated with a classically for- frequently employed in one-dimensional models of driven 

bidden tunneling process. It is, on the other hand, also atoms (e.g., |38|), may not be permitted, since it allows 

clear that a full-blown quantum treatment of such highly the electrons to penetrate to the other side of the nucleus, 

correlated states under external driving is still beyond what inevitably destroys the frozen-planet configuration, 
reach of the presently available computational facilities. 



For unperturbed helium, where the total angular momen- 
tum is conserved, autoionizing states can presently be cal- 
culated for quantum numbers of the inner electron reach- 
ing N ~ 15 ... 20 Therefore, in the near future 
there is not much hope for an accurate treatment of the 
periodically driven problem, with a large number of angu- 



Consequently, we write the Hamiltonian that generates 
the quantum dynamics of the driven collinear configura- 
tion as 



H 



1 d 2 


i d 2 


Z 


Z 


1 


2~d4 ~ 




Z\ 




H 

z\ - z 2 


F . 

— smut 

LU 


n d 

\i dzi 


1 d N 
i dzi 


)- 


F st (zi+%),(4) 
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where the exact Coulomb interactions between the charged what will turn out as a convenient form for our further 
particles are taken into account. Here, Z\ and z 2 represent treatment, 
the Cartesian (not spherical) coordinates of the electrons 

along the field polarization axis, with < z 2 < Z\ < oo 3.2 Floquet theory and complex scaling 

- i.e., both electrons are located on the same side of the 

Due to the temporal periodicity of the Hamiltonian, the 

nucleus, with electron 2 closer to it than electron 1. Note 

Schrodinger problem represented by (|6|) is conveniently 

that the electrons effectively appear as distinguishable par- 
treated in the framework of Floquet theory Any 

tides. This ansatz is consistent with the restriction to 

solution of the Schrodinger equation can be expanded in 

collinear motion, since in the classical collinear configu- 

a set of 27r/w-periodic quasienergy-eigenfunctions ?/> t , 

ration the electrons are not able to pass each other and to 

exchange their positions. Furthermore, it represents quite ^* = ^ d£Cei\)\ ] e 8 , with V't+Lr/w = W eC 
good an approximation also for the three-dimensional frozen- 

(7) 

planet configuration, since the latter is characterized by a which satisf y the Floquet eigenvalue equation 

very small overlap between the outer and the inner elec- ( TT . 9\ is) «?/(£) / D \ 

I H-i— J ipl = £ ip\ i (8) 

tron's densities, and therefore by a very small splitting 
between singlet and triplet states Jn| . Precisely as for the 

oo 

classical Hamiltonian, the external electric field consists ip( £ ' = l ip^e lku ' t , (9) 

of an oscillating part F - which, in contrast to ([!]), is in- 
corporated in the velocity gauge in order to ensure better 
convergence of the numerical calculation j3i| - and of a 
static component F st > which forces the electrons away 



A Fourier series expansion of tp\ , 



k— — oo 

yields a time-independent eigenvalue problem for the Fourier 
components ip k > where the effective Hamiltonian matrix 



from the nucleus. 

In a last step, we introduce the coordinates 



contains the stationary part of the Hamiltonian shifted by 
kto as diagonal, and the Fourier components of the peri- 
odic driving as off-diagonal (block) elements. This results 
in an w-periodic spectrum of quasienergies £. 
£ = z i ~ z 2 , V = Z2 , (5) p or atomic systems, the Floquet spectrum is absolutely 

continuous: each bound state of the unperturbed atom is 
coupled to the atomic continuum via multiphoton transi- 
tions, and therefore appears as a resonance structure in 
the spectrum - in analogy to autoionizing states in mul- 



which independently lie in the range < £,77 < 00. In 
these perimetric coordinates [ [f0| , the Hamiltonian is rewrit- 
ten as 

d 2 1 d d 2 Z Z 1 



77 rj £, tiply excited atoms, which are coupled to the continuum 

F Id 

- — ■ sin tot y— - Fst (£ + 277), (6) via configuration interaction Q (see, e.g., Q for a case 
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in triply excited lithium) . In order to separate these res- composed of the real- valued Sturmian functions 
onances from the flat background of the continuous spec- 



trum, we use the method of complex scaling [^5[^,^|. ' ^ \fn a \ a 



sfHo = t^ 2 A exp {L) L M 1 (*A ] , (14 , 

\/n a \ a J \ a 



It consists in the complexification of coordinates and mo- 
menta (not of time!) according to [] 



where the denote the associated Laguerre polynomi- 
als p9|. The real parameter a > fixes the length scale 



ig . d . d _ ig . ( |l3| ) exhibits two independent scaling parameters, a and 

/3, which permit the optimal tuning of the basis with re- 



d d 

£ y £ e ie —i > _j e ~ ie of the Sturmians. As a consequence, the product basis 

,d_ _.d_ 

dr/ dr/ 

through application of the nonunitary complex scaling op- 

spect to the relative excitation of the electrons. Also note 

erator R(9) on ibi £ \ . r a \ , . 

that the Sturmians 6n scale at least linearly with £, for 

R(9)ip[ £ \^ rf) — e ie tp[ £ \^e ie r]e ie ) =■ ip^ £ \^ rf) (11) £ 0- This implies that the Coulomb singularities of the 



Accordingly, the Hamiltonian is transformed as 



Hamiltonian (||) do not lead to divergent matrix elements 
in the product basis jl^). 

H — > H e = R(6)HR{-6), (12) Multiplication of (§) by £77 leaves + 77) as the only 

non-polynomial contribution in £ and 77, after inserting 
(§). We factor out (£ + 77) fr om the eigenvectors of (g), 



# } (^) =: (£ + r7)^ £) (£,77), (15) 



what turns (g) into a complex symmetric rather than her- 
mitian eigenvalue problem, with eigenvalues in the lower 
half of the complex plane. In particular, resonances of the 
"real" eigenvalue problem (i.e., at 9 = 0) appear, for finite 

9, as isolated, discrete complex eigenvalues 8 = E — ir/2: and expand <jyf rf) in the Sturmian product basis. This 
their real and imaginary parts correspond to the energies factorization is justified since the triple collision between 
E and to the half widths at half maximum r/2 (HWHM) the electrons and the nucleus - i.e., the event that cor- 
of the resonances, respectively. responds to £ = 77 = - is classically suppressed in the 

collinear frozen-planet dynamics (in contrast to the eZe 
3.3 Expansion in Sturmian basis functions configuration with the electrons on opposite sides of the 

nucleus); hence, the quantum probability near £ = 77 = 



The complex-scaled Floquet Hamiltonian is now expanded 
in the product basis 



is very small. Note that, due to normalization, (15) entails 
a different transformation law for <f> t (£,77) as compared 
{S^i^S^irj) :n,m>l} (13) to # } (£, 77), under the action of R{9) 



1 Note that, due to (^), complexification of £, 77 and of Zi, Z2 
is equivalent. R(9)^ £) (£, 77) = e 2l9 # } (£e W , V e* e ) =: ^(£,77) . (16) 
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Altogether, the generalized Floquet eigenvalue equa- the Lanczos algorithm with inverse iteration [|50|. Conver- 
tions finally reads gence of the complex eigenvalues is contolled by increasing 



><ZvVe(Z + v)i$,L 1 (Z,v)-$,l +1 (Z,v)) = 0, (17) 



*(£) 

the number of Fourier components 4> g 1 ana - °f the basis 



functions (£) , Sm (v) 1° be retained in the diagonal- 
ization, as well as by variation of the complex scaling angle 
with the complex scaled operators H oe , V e given by 9, and of the length scale parameters a and (3 which spec- 

/ g I q q2 \ ify the (truncated) basis set ( jfjj ) used in the calculation. 

All numerical results presented hereafter are converged to 



"<9£ 2 2dr] 2+ d£dr)J 

Z Z V 

— + -r I e 



£ + 77 rj £J machine precision. 

f,iK + 2l)e", (IS) 



and the 4> s £ l(S,, rj) the Fourier components (see (||)) of the 

(£) 

wavefunction tp g t (£,77). Introducing the scalar product 

oo />oo 



3.4 Visualization of the wavefunction 



The energies E and widths r extracted from the com- 



' I I £77 f(C r l)9(^i r l)d^dr], (20) eigenvalues £ allow to classify the Floquet states in 



Jo 



.,, , , v. v .v rij. j t 1 series and to determine their stability; in general, how- 

with respect to which the bturmians are orthogonal, we J ' 

1 , . , . , 1 , ever, they do not provide any direct information about 

obtain a complex symmetric eigenvalue problem. ' J r J 

. , , i j 1 , • ,1 j the localization properties of these states in configuration 

All the opertor-valued terms acting on the 1 and r ^ ° 

l(£) (FT^K u j 1 . , , ,, or in phase space, which in themselves can carry impor- 

<pg tj-i m (17) can now be expressed as polynomials ol the 

, 1 j j ; r ,1 0i i- /FT7i\ tant physical information. To extract the latter from the 

elementary ladder operators of the Sturmian basis (14J) r J 

/ icra ittj > TT t ,, „ . eigenstates of the complex symmetric eigenvalue problem 

(see [|23|,|49|). Hence, upon expansion of the Courier com- r j o t 

1(E) ■ ai. Di i i i_ jtt]\ (171), some care has to be taken, due to the non-unitarity 

ponents <p e j. m the bturmian product basis (|l_^j , K > — v ' ' J 

of the dilation operator R(6). The technical and compu- 



&V) = E ci n%StH0S^(v), (21) 



tational details of this procedure are described in detail in 

we obtain strict selection rules which permit nonzero cou- ©HI > and we onl y im P ort here those results of immedi- 

plings between the coefficients C&% and C^, only if ate relevance for our specific purpose. 

\n — n'\ < 4, I to — to'| < 4, and \k — k'\ < 1. The eigen- We first need the density of the physical wave function 

value problem thus becomes a sparse, complex symmetric (at real energy E) ipl (zi,Z2), which is obtained (up to 

matrix equation with rather narrow bandwidth, what al- a normalization constant), for a well-isolated resonance 

lows to apply efficient diagonalization methods based on energy £ = E—ir/2 (i.e., \£j— £\ ^> f ,V£j 7^ £), from the 
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image of tpg j under the inverse dilation operator R(—9), of the driving field), and at fixed position z% ~ of the 



inner electron. The squeezing parameter oj s in this defini- 
tion determines the resolution of the Husimi distribution 
in position and momentum, respectively. Large values of 
Lo s yield good resolution in position space and bad resolu- 
tion in momentum space, and vice versa for small values of 
\tJj\ £ \zi, Z2)\ 2 — Re(zi, Z2\R{— 9)\^g £ ^)(zi, Z2\R{— 0)\ijjg £ ^ t ) . to s . For quantum states associated with the frozen-planet 

(23) configuration, the most appropriate choice for lo s is given 
Correspondingly, the Husimi representation Q(z 1 , pi) ||| by the intrinsic frequency scale wj of the configuration 



R(-8)4 £ J = R(-6) z 1 ^l(z 1 -z 2 ,z 2 )e ik ^ t , 

k— — oo 

(22) 

together with the image of the corresponding left eigen- 
vector 4 £) _ t [HIM 



|5l| in the classical phase space coordinates of the outer (see Section 2.1) - this is the value at which the Gaussian 



electron (subscript 1) is obtained by projection on Gaus- wave packet Xp(l) optimally fits the ground state of the 
sian wave packets Xpi( z i) which are centered at position effective potential experienced by the outer electron. We 
z\ and propagate with average momentum py. therefore choose u> s = 0.3 N^ 3 (see Fig. ||), with N the 

inner electron's principal quantum number of the state to 



Q{z x ,p x ) ^{ Xp M)\R{-^f}){x P A^)\R{-mfl t , , 

(24) 

where the bar denotes complex conjugation. 

In order to establish a precise analogy with our above 
construction of the two-step Poincare section of the classi- 



be visualized (note that TV remains approximately a good 
quantum number in the presence of the external field, if 
the latter is not too strong). 



cal dynamcis (see Section 2/1), we define the overlap ma- 
trix element in (E4|) as 



4 The unperturbed collinear atom 



(XM\R(-0)\4 £ 1)= dz[{z[,4 o) \R{-0)\^) 

Jo With the above theoretical machinery for our quantum 

x exp { ^uj s (z 1 z x ) ipiz^J ■ (25) treatment of the driven, collinear frozen planet configu- 

. . . . . ; . . ration at hand, we can now study in detail its physical 

Ihus, the quantum phase space probability density is eval- 

, , .. nii , properties. As a "warm-up", and equally so as a first as- 

uated at fixed time t (corresponding to a fixed phase uit 



2 Only the time-dependent part of the wave function is com- sessment of the physical implications of the restricted di- 

plex conjugated to transform right into left eigenvectors, since mensionality of our problem, we start out with the spectral 

the Fourier components %j)fl of the eigenvectors obey the nor- properties of the unperturbed (F — F st — 0) frozen planet 

malization condition for complex symmetric rather than her- configuration confined to a single dimension of configura- 

mitian matrices [§,||]! tion space. 
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Fig. 9. Energies and ionization rates JT/2 (HWHM) of the autoionizing states of collinear helium belonging to the series labeled 
by the inner electron's quantum number TV = 8. Note the relatively small widths of the lowest members of the series, which we 
attribute to the fact that these states are fully localized within a regular, bound domain in phase space. The coupling mediated 
by a driving field with frequency uj = 0.2 iV -3 is sketched by the arrows of length uj. The particularly strong, near-resonant 
coupling between the second and the third state (at E — —0.034 and E = —0.0336, respectively) gives rise to the nondispersive 
two-electron wave packet associated with the classical 1:1 resonance (see Fig. |l|). 

4.1 Energies and phase space distributions ber arises from the quantization of the outer electron's 

oscillation within the effective potential, which, as in one- 



electron atoms, gives rise to Rydberg series and continua, 
due to its asymptotic \ jr dependence. 



Corresponding to the fact that its classical counterpart 
is almost completely integrable (Fig. f|(a)), the quantum 
spectrum of the unperturbed collinear frozen-planet con- 
figuration is characterized by typical signatures of regular 

dynamics. In particular, each autoionizing state can be As an example, we show the energies and widths of 
classified by two quantum numbers which are associated the autoionizing states of the series labeled by N — 8 in 
with the two separating modes of the classical dynamics, Fig. ^. We clearly recognize a regular Rydberg progres- 



as discussed in Section 2.1 (in contrast to the chaotic eZe sion towards the N — 8 threshold at E = —0.3125 a.u. 

configuration |32| , where such a classification breaks down States belonging to different series are energetically well 

already for moderate double excitations). The quantiza- separated: The lowest quantum state with A = 9 lies at 

tion of the inner electron's Kepler mode gives rise to the E ~ —0.027 a.u., and the Rydberg series of AT = 7 termi- 

quantum number N > 1. For each N, we obtain a series of nates at E ~ —0.04 a.u. As was verified numerically, an 

autoionizing states (or bound states, for N = 1) which are energetic overlap of series belonging to different quantum 

labelled by the quantum number n. This quantum num- numbers A^ occurs only above N = 20. 
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Fig. 10. Wavefunctions of the first three eigenstates of the N = 8 series: (al, bl) lowest (n = 1), (a2, b2) second (n = 2), (a3, 
b3) third state (n = 3) of the series, (al - a3) Probability densities as a function of the coordinates of the outer (zi) and of 
the inner (22) electron; (bl - b3) Husimi densities as a function of position and momentum of the outer electron (the contour 
lines are plotted on a linear scale). In the coordinate of the inner electron, the wavefunctions are essentially given by the JVth 
hydrogenic eigenstate of the unperturbed, one-dimensional atom, whereas in Z\ they display different excitations within the 
effective potential (see Fig. |^). Their Husimi distributions clearly show that they are well localized along regular tori in the 
underlying classical phase space (see Fig. ^a)). 



The near-separability of the two modes associated with 
the electrons is furthermore reflected by the probability 
density ( |23| ) of the eigenstates (with k = in (p2]), due 
to the time- independence of the problem). In the coordi- 
nate zi of the inner electron, the wavefunctions are essen- 
tially given by the ./Vth hydrogenic eigenstate of the one- 
dimensional atom (most reminiscent of extremal parabolic 
states of three-dimensional atomic hydrogen), whereas in 
z\ they display distinct excitations within the effective po- 
tential governing the outer electron's motion. In Fig. [l0|, 
the fundamental (al), the first excited (a2), and the sec- 
ond excited state (a3) of the effective potential (for N — 
8), labeled by the outer electron's quantum numbers n = 
1, n = 2, and n = 3, respectively, are plotted in configu- 
ration space. The intimate correspondence of these eigen- 



states with regular structures of the underlying classical 
phase space is highlighted by their Husimi distribution 
( pi[ ) shown in Fig. |lO](bl - b3). Comparing these phase 
space projections of the quantum eigenfunctions with the 
corresponding Poincare surface of section (Fig. |](a)) un- 
ambiguously illustrates the neat localization of the eigen- 
functions along invariant tori of the classical dynamics. 

4.2 Ionization rates 

The "regular" character of the wavefunctions is also man- 
ifest in the ionization rates r/2 (HWHM) of the associ- 
ated spectral resonances. We see in Fig. ^ that for the 
N = 8 series these widths first increase with n, and then 
decrease again towards the ionization threshold. While 
this latter decrease is attributed to general scaling laws 
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Fig. 11. Energies (a) and ionization widths (c) of the lowest autoionizing eigenstate n — 1 (with the outer electron localized 
at the minimum of the effective potential of Fig. ^]), plotted as a function of the inner electron's quantum number N. We see 
that the ionization rates decrease near-exponentially for N > 4, what suggests that the associated eigenstates are fully localized 
in the bound part of phase space, and only decay via tunneling through the classical phase space barriers that protect the 
frozen-planet configuration against classical autoionization. These rates of the collinear model are compared to the energies (b) 
and decay rates (d) of the analogous frozen-planet eigenstates in three-dimensional helium, published in Ref. [[L9|], for singlet (+) 
and triplet (O) symmetry. While the energies agree quite well with their collinear counterparts (a, b) , the ionization rates of the 
three-dimensional frozen-planet states (d) lie by several orders of magnitude above the widths of the corresponding eigenstates of 
the one-dimensional model (b) (note the different scales of the vertical axes in (b) and (d)). This indicates that the autoionization 
of the three-dimensional frozen-planet states is mediated by the transverse degrees of freedom of the configuration. 

of r in the Rydberg regime of atomic one-electron excita- a function of the inner electron's quantum number N . In 

tions HI , the extraordinary stability of low n states is a Fig. |ll|(c) , we observe that the width of this lowest eigen- 

consequence of their full localization in the bound part of state in the effective potential of Fig. |^ first increases for 

phase space. A coupling to the unbound, ionizing domain N < 4 (due to imperfect localization of the eigenstate 

is then only possible via a classically forbidden process - on classical phase space structures in the deep quantum 

tunneling through the phase space barriers that confine regime) and then decreases exponentially with N. This 

the frozen-planet configuration. exponential law represents a clear signature of quantum 

tunneling: The rate associated with such a process gener- 

This latter mechanism is beautifully illustrated in a „ , , • n .. -, , -, . , . , 

J ally decreases exponentially with the semiclassical action 

plot of the ionization rate of the lowest state (n — 1) as 
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Sn across the tunneling barrier, and this latter quantity On the other hand, the energies of the autoionizing 

scales, as any other action variable, linearly with N, due states of three-dimensional helium agree quite well with 

to the general scaling laws (||) - i.e., Sn = NSi, with their collinear counterparts, as apparent from Figs. p"T|(a) 

Si ~ 0.8 directly extracted from Fig. [ll](c). and [ll](b) (apart from a small shift due to the additional 

For comparison, we plot in Fig. |ll](d) the ionization quantization of the transverse modes in the 3D configura- 

rates of the near-collinear frozen-planet states in three- tion jL9|). This supports our expectation that the essen- 

dimensional helium (with minimal excitations in the trans- tial structure of the eigenstates of the three-dimensional 

verse degrees of freedom) , as reported in . Also these frozen-planet configuration is well reproduced by the collinear 

widths decrease on average exponentially in N, with a model, 
scaled action S± not too different from our one-dimensional 

model (fluctuations with respect to this exponential de- _ _. ... 

5 The collinear atom under periodic driving 

crease arise from the mixed regular-chaotic phase space 

structure of three-dimensional helium, see also the discus- 5.1 Nondispersive wave packets in the Floquet 



sion in Section 5.2). However, the three-dimensional rates spectrum 



are systematically enhanced by several orders of mag- 

For our investigation of the collinear frozen planet con- 

nitude as compared to the one-dimensional model! This 

figuration under external driving, we employ the field pa- 
suggests that the three-dimensional helium dynamics pro- 

rameters w = 0.2 N~ 3 and F = 0.005 N~ A which induce 

vides efficient decay channels associated with the trans- 

the classical phase space structure shown in Fig. ^|, where 

verse degrees of freedom, which are not incorporated in the 

both the intrinsic island and the 1:1 resonance are very 

restricted one-dimensional dynamics of our model Hamil- 

well pronounced. The coupling induced by this driving 

tonian (|).| 

— field is illustrated by the arrows in Fig. H for the N = 8 

As a side remark, this observation also bears a caveat as 1-1 

, , . „ , , . . . series. We see that a particularly strong interaction is in- 

oversimplmed models of the three body Coulomb problem are 

j £, ,-, , , . i ... , in , , duced between the second (n = 2) and the third (n = 3) 

concerned: often the Coulomb singularities in such ID model v ' v ' 

systems are smoothed such as to prevent the electrons (and state > since the energetic spacing between these states is 

possibly the nucleus) from head-on collisions §jj,|§, which, very close to the driving frequency. Provided the phase 

on a first glance, could be most desastrous for the stability space volume of the 1:1 resonance island is large enough to 

of the ID dynamics and, hence, enhance its decay probability, support a fully localized quantum state (which is indeed 

The present comparison of exact (no approximations beyond the case for N — 8, as can be semiclassically estimated 



dimensionality) ID and 3D quantum calculations shows that namical stabilization effects which result in such impressively 
this is too simplistic an argument and completely neglects dy- stable quantum objects as the frozen planet. 
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by the EBK quantization criterion), this near-resonant Comparison with Figs. [lO|(al - a3) shows that the eigen- 
one-photon coupling gives rise to the nondispersive wave function is essentially composed of a time-periodic super- 
packet associated with the classical nonlinear resonance. position of the low-lying eigenstates of the N — 8 series. 

tt FT71/ \ i i, i • • While nearly stationary in the coordinate of the in- 

b lg. |12|(a) shows the evolution of the quasienergies as- J 

. , , .., n , tt, , . r , , ner electron, the wave packet oscillates back and forth in 
sociated with the hrst three Moquet eigenstates of the 

tit q . r r r. ii ,., , „ „ zu between the turning points of the effective potential 

N = 8 series, as a lunction of the held amplitude t (for ' 01 r 

fixed driving frequency lu — 0.2 A~ 3 , and F st = 0). In the (^& - @)' 

limit .F = 0, the energies E correspond to the unperturbed Quite naturally, for the same field amplitude and fre- 
autoionizing levels folded into the Floquet zone due to the quency as above, the diabatic continuation of the low- 
w-periodicity of the Floquet spectrum. Switching on the est state n = 1 of the N = 8 series represents the Flo- 
driving field causes level shifts due to the repulsion be- quet state that is anchored to the remainder of the un- 
tween strongly coupled, near-resonant Floquet states. A perturbed regular structure of the frozen-planet config- 
particularly significant case is encountered for the state uration. Fig. fl4|(al - a3) shows its Husimi distribution 
n = 2, near-resonantly coupled to the states n — 3 and which is fully localized on the intrinsic island and re- 
ft = 1. Its quasienergy exhibits a pronounced shift at al- mains nearly stationary over one field period. Further- 
most constant, positive slope. This already indicates the more, Floquet states corresponding to higher quantum 
wave-packet character of the associated eigenfunction in numbers n > 2 of the unperturbed system are entirely 
configuration space, since it reflects a large dynamic dipole associated with the chaotic phase space domain. As an 



moment of the electronic density (proportional to the ex- example, Fig. 14 (bl - b3) shows the Husimi density of 

pectation value of the electronic dipole z\, averaged over the diabatic continuation of the state n = 3, which is 

one cycle of the drive). dominantly localized on the unstable fixed point of the 

rp, iii , r ,,. , ii • lj.* 1:1 resonance (compare Figs. [14](b3) and (c3)), and along 

Ihe wave packet character of this latter state is ulti- v ^ ° <— > y ' \ i n t> 

, i n , , .. i • i • i i the associated separatrix layer, 
mately confirmed by its Husimi density which is plotted 

in Fig. [l3](bl - b3) for F = 0.005 TV -4 , at driving phases A similar scenario is encountered for different quantum 

u)t = 0, 7r/2, and tt. We see that the Floquet eigenfunc- numbers N of the inner electron (with adjusted field pa- 

tion is fully localized on the 1:1 resonance island in phase rameters ui = 0.2 N~ 3 and F — 0.005 A^ -4 ). Starting from 

space, and that it precisely follows the classical time evo- A^ = 5, we find, for each N, one Floquet state that is lo- 

lution. The wave packet dynamics in configuration space calized on the intrinsic island (the diabatic continuation of 

is visualized in Fig. [l3|(al - a3) where we plot the proba- the island's ground state in the unperturbed spectrum), 

bility density of the state n — 2 as a function of z\ and z%. as well as one nondispersive wave packet state localized 
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Fig. 12. (a) Evolution of the first three states of the N — 8 series in the Floquet spectrum, for driving frequency uj = 0.2 TV~ 3 . 
The quasienergies E are plotted as a function of the field amplitude F. The diabatic continuation of the second state (n = 2, 
upmost circles) undergoes a particularly pronounced, almost linear shift with F. This is consistent with the wave packet character 
of the Floquet state into which the unperturbed n = 2 state evolves (see Fig. ^). (b - d) Ionization rates _T of the diabatic 
continuations of the n — 1 (d), n = 2 (b), and n = 3 state (c). We see that the rates globally increase with increasing F. Local, 
pronounced maxima occur at avoided crossings between these states and energetically higher-lying eigenstates of the N = 8 
series in the Floquet spectrum, as apparent from a comparison with (a) (with the diabatic continuations of n — 1 and n = 3 
marked by the encircled eigenvalues in the lower and the middle part of (a), respectively). 



on the 1:1 resonance island, represented by the diabatic 
continuation of the second unperturbed state for N < 10, 
and by the continuation of the third or a higher state for 
N > 10. As for N = 8, these wave packet states are char- 
acterized by a pronounced, almost linear level shift with 
increasing F towards higher quasienergies (Fig. |l5|) . 

Fig. |l^ compares the Husimi distributions of the wave 
packet states for N = 5 and N — 15. Since the phase space 
area of the island increases linearly with TV while h remains 
constant, the electronic density is, with increasing N, more 
and more localized around the center of the 1:1 resonance 



island. At N — 15, the contour lines of the Husimi density 
reproduce the classical island structure in the Poincare 
section already very well. In fact, at such high values of 
TV, the island is already large enough to trap an additional 
wave packet state, giving rise to a first excited mode of 



the nondispersive wave packet (see Fig. 16 (cl - c3), and 
Ref. for an analogous case in driven hydrogen). 

5.2 Lifetimes of the wave packet states 

As already mentioned in the Introduction, the time scale 
over which the nondispersive wave packet follows the res- 
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Fig. 13. Nondispersive two-electron wave packet of the N = 8 series, for uj = 0.2 N~ 3 and F = 0.005 N . We plot the 
probability density in configuration space (al - a3), and the Husimi density in phase space (bl - b3), of the diabatic continuation 
of the state n = 2 (see Fig. [li]), at driving field phases (al - cl) uit — 0, (a2 - c2) ty/2, (a3 - c3) n. The corresponding Poincare 
sections of the classical phase space are plotted in (cl - c3). We see that the Floquet wavefunction is well localized on the 1:1 
resonance island, and faithfully tracks its classical time evolution. 

onantly driven classical orbit without spreading is only neling barrier that confines the resonance island at N = 1. 
limited by the finite ionization rate of the corresponding In contrast to the unperturbed collinear configuration dis- 
Floquet state. Since this ionization proceeds via a clas- cussed in Section || (Fig. |ll](c)), however, such an expo- 
sically forbidden tunneling process, the associated decay nential decrease is not observed in our calculations. As we 
rate is, as in driven hydrogen |]6|,p[,[Tc|], very small. For the see in Fig. [l7| the ionization rates exhibit, in the range 
wave packet state of the N — 8 series {uj = 0.2 7V -3 and N = 5 ... 15, strong, quasi-erratic fluctuations between 
F = 0.005 N~ 4 ) for instance, we obtain Fj 2 ~ 1.77- lO" 11 T/2 - 10" 12 a.u. and r/2 - 10" 8 a.u., corresponding to 
a.u., which corresponds to a lifetime of r ~ 1.8 • 10 6 field lifetimes between 10 4 and 10 7 field cycles. No significant 
cycles before the wave packet population is appreciably tendency towards lower values for increasing N can be 
redistributed over the ionization continuum. unambiguously identified. 



In general, one expects the ionization rate of the wave 
packets to decrease on average exponentially ~ exp(— SN) 
with TV, where S is the semiclassical action across the tun- 



The observed fluctuations of r represent a clear signa- 
ture of the mixed regular-chaotic structure of the classi- 
cal phase space (see, e.g., Fig. |l6|(dl-d3)). In contrast to 
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Fig. 14. Husimi distribution of the states n = 1 (al - a3) and n = 3 (bl - b3) of the N = 8 series, for lo = 0.2 TV -3 and 
F = 0.005 iV -4 , at driving phases wi = (al,bl), 7r/2 (a2,c2), it (a3,b3), compared to the corresponding classical Poincare 
sections (cl - c3). We see that the diabatic continuation of the lowest state of the N = 8 series is fully localized on the intrinsic 
island and remains essentially stationary over one field cycle. The diabatic continuation of the third state, on the other hand, is 
associated with the chaotic phase space domain. More precisely, it is anchored to the unstable fixed point of the 1:1 resonance 
and stretches out along the associated separatrix structure (see also pLlM for analogous states in driven hydrogen atoms). 



globally regular systems such as the unperturbed Zee con- 
figuration (see Section [l^), the coupling from the locally 
bound domain to the ionization continuum is not given by 
a single tunneling rate. Instead, it is mediated by a two- 
step process which, besides tunneling through the barri- 
ers of the regular island, involves also chaotic transport 
through irregular phase space outside the island towards 
the continuum. It is in fact this latter, chaotic component 
of the coupling process that generates the fluctuations of 
the decay rates, on top of their on average exponential 
decrease with N . This chaos- assisted tunneling mecha- 



nism, originally formulated for tunneling transitions be- 
tween symmetric regular islands in bounded systems | p6j , 
was indeed proven to govern the decay of nondispersive 
wave packets in driven one-electron atoms fl2|| , the lat- 
ter exhibiting similar quasi-erratic deviations of r from a 
purely exponential law (see also jp7fl). 

In the Floquet spectrum, the chaos-induced modifi- 
cation of the tunneling process can essentially be under- 
stood via the influence of near-resonant chaotic states on 
the ionization rate of the wave packet state. As we see 
in Fig. [T^(b-d), the widths of the energetically low- lying 
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Fig. 15. Energies and widths of the nondispersive wave packet at u) = 0.2 TV -3 for TV = 9 (al, a2) and TV = 10 (bl, b2), as 
a function of the scaled field amplitude F/N~ 4 . (al, bl) Evolution of the quasienergies E in the Floquet spectrum. The 1:1 
resonance state (marked by circles) diabatically originates from the unperturbed state n = 2 at TV = 9, and from the state 
n = 3 at TV = 10. As for TV = 8 (see Fig. |l^a), its quasienergy undergoes a significant shift towards higher values of E with 
almost constant slope, what indicates the pronounced wave packet character of this eigenstate in configuration space. (a2, b2) 
Ionization widths of the wave packet states. As in Fig. |T^(c) , local, sharp enhancements of F occur near avoided crossings 
with higher-lying states of the same series, which are associated with the chaotic phase space domain. Notice that the precise 
positions and strengths of the avoided crossings and, consequently, their influence on the wave packet's lifetime at a given scaled 
field amplitude F/N~ 4 , can be completely different for TV (al, a2) and TV + 1 (bl, b2). 



members of the TV = 8 series generally increase with in- 
creasing F. Local, pronounced enhancements of 7T 1 , how- 
ever, occur near avoided crossings (corresponding to mul- 
tiphoton resonances) with higher-lying states of the same 
series, which are associated with the chaotic phase space 
domain (compare Fig. |l2|(b-d) with Fig. [l2](a)). Hence, 
at a given scaled field amplitude, e.g. F = 0.005TV~ 4 , the 
decay rate of the wave packet state depends sensitively 
on the effective interaction with nearby chaotic states in 
the spectrum, the latter being essentially determined by 
the energetic distance to those states, as well as by the 
corresponding coupling matrix elements. Since these in- 



gredients vary in a random-like way with the quantum 
number TV of the inner electron (at fixed F/N~ 4 ), the 
resulting effect on the wave packet's ionization rate may 
turn completely different when changing TV to TV + 1 (com- 
pare here Fig. |l^(a,b) to Figs. |l5](al,a2) and (bl,b2)), and 
therefore cannot be predicted by simple semiclassical con- 
siderations. 

The average exponential decrease of the ionization rate, 
on the other hand, should in principle be amenable to 
semiclassical estimations, as it is essentially determined 
by the height of the tunneling barrier that confines the 
resonance island in phase space. Semiclassical approaches 
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Fig. 16. Nondispersive two-electron wave packet for N = 5 (al - a3) and N = 15 (bl - b3), as well as the first excited mode 
of the nondispersive wave packet for N = 15 (cl - c3), at U3 = 0.2 iV~ 3 and F = 0.005 iV -4 . The Husimi density of the Floquet 
states is plotted for the phases (al - dl) cot — 0, (a2 - d2) tt/2, (a3 - d3) n, and compared to the corresponding classical Poincare 
sections (dl - d3). With increasing N, the wave packet is more and more localized on the center of the 1:1 resonance, due to 
the decreasing relative size of 2nh (as compared to the typical atomic scales, see Eqs. (Ej)) indicated by the dashed rectangles in 
(al) and (bl). At N = 15, the contour lines of the Husimi density reproduce the structure of the island in the classical Poincare 
section rather well. 

to quantitatively evaluate the associated tunneling rates 

for nonintegrable systems are presently under construe- 

well for N = 8,9, 10, as equally observable in Figs, [l^ and |l5| 

tion (see ]58fl).p| 

This is a signature of the transition from globally regular to 
4 Also note, that the critical field amplitude needed to in- dominantly chaotic dynamics in the classical phase space, and 
duce an appreciable increase of the decay rates beyond their highlights the relevance of the classical scaling laws (||), also 
minimum value at small field amplitudes coincides remarkably for the quantum mechanical ionization process. 
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Fig. 17. Ionization rates (HWHM) of the wave packet state localized on the 1:1 resonance island, as a function of the inner 
electron's quantum number N (uj = 0.2 TV -3 and F = 0.005 N~ 4 ) , for different strengths F s t of the static field. Circles connected 
by thick solid lines: F Bt = 0; squares connected by thin solid lines: F st = 0.001 iV -4 ; diamonds connected by dotted lines: 
F Bt = 0.002 7V~ 4 . The strong fluctuations of the rates are attributed to the coupling to Floquet states that are associated with 
the chaotic domain of phase space. 

5.3 Influence of an additional, static electric field duction of their size, the two dominant regular islands are 

not significantly affected by the static perturbation. 

Let us finally discuss to which extent the wave packet 

dynamics of the collinear atom is modified by the applica- Accordingly, we find that the phase space projections 
tion of an additional, static electric field with a strength of the Floquet states localized on these islands remain es- 
of up to 20 percent of the driving field amplitude. As al- sentially unmodified, too, in presence of the static field. 



ready pointed out in Section 2.3, such a static field is This is illustrated in Fig. 18(bl - b3) which shows the 



required to stabilize the driving-induced nonlinear reso- resulting Husimi density of the nondispersive wave packet 

nances in the classical dynamics of the three-dimensional state of the N = 8 series at w — 0.2 N~ 3 and F = 

atom. While the transverse dynamics of the driven con- 0.005 A^~ 4 . However, while no visible change occurs in the 

figuration is therefore substantially modified by the static phase space localization properties, the ionization rate is 

field, the structure of the collinear phase space remains significantly affected by the static field: it increases from 

nearly unaffected. This is seen in Fig. [jj(al - a3) which T/2 ~ 1.77 • lO" 11 a.u. at F st = 0, to r/2 ~ 3.05 • I0~ 9 

shows Poincare sections of the collinear phase space at a.u. at F st — 0.002 A^~ 4 , corresponding to a decrease of the 

lu = 0.2 and F = 0.005, for gradually increasing static wave packet's lifetime from 1.8 • 10 6 to 1.0 ■ 10 4 field cycles, 

field components F st — . . . 0.002. Except for a slight re- Since the static field notably affects the asymptotic part 
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Fig. 18. (al - a3) Poincare section of the classical dynamics at u = 0.2 and F = 0.005, with a static field Fst = (al), 
0.001 (a2), 0.002 (a3). (bl - b3) Husimi density of the wave packet state localized on the 1:1 resonance island for N = 8, at 
u = 0.2 iV~ 3 and F = 0.005 N~ 4 , calculated in presence of a static field F at = (bl), 0.001 N~ 4 (b2), 0.002 iV' 4 (b3) (at 
driving phase uot = 0). The phase space projection of the wave packet remains almost invariant when increasing the static field 
from Fst = to 0.002 iV -4 , in agreement with the observation that also the classical 1:1 resonance island is largely unaffected 
by the additional perturbation. 

of the effective potential experienced by the outer electron linear resonances, which actually reflect a phase locking 
(see Fig. ^|), and thereby modifies the transport towards phenomenon |l6| ] between the unperturbed quantum dy- 
the continuum, such an enhancement of r with F st is not namics and the driving field, allow to launch non-dispersive 
unexpected. However, Fig. [l^ shows that the ionization two-electron wave packets along the highly correlated frozen 
rates do not always increase with increasing F st . Instead, planet orbit of collincar helium, with tiny, tunneling in- 
a rather erratic dependence of r on F st and N is encoun- duced leakage towards the atomic continuum. As already 
tered, what calls for a more quantitative (and possibly observed in simpler one electron systems under external 
statistical) analysis of the problem. driving, these wave packets can be identified by the char- 

acteristic parameter dependence of their energies and de- 
cay rates, with strong signatures of chaos assisted ioniza- 
tion. Our ID quantum calculations corroborate an inti- 
mate quantum-classical correspondence, already at rather 
weak excitations (N > 5) of the inner electron, even in the 
presence of a static electric field component added to the 



6 Conclusion 

Our present study provides abundant evidence for the sta- 
bilizing effect of classical nonlinear resonances on period- 
ically driven quantum systems, in the presence of strong 
electron-electron interactions. We have seen that such non- 
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